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Abstract 

The area of the unitarity triangle is a measure of CP-violation. We introduce the 
leptonic unitarity triangles and study their properties. We consider the possibility of re- 
constructing the unitarity triangle in future oscillation and non-oscillation experiments. 
A set of measurements is suggested which will, in principle, allow us to measure all 
sides of the triangle, and consequently to establish CP-violation. For different values of 
the CP-violating phase, 5d, the required accuracy of measurements is estimated. The 
key elements of the method include determination of jC/esI and studies of the — 
survival probability in oscillations driven by the solar mass splitting Am-g^^. We sug- 
gest additional astrophysical measurements which may help to reconstruct the triangle. 
The method of the unitarity triangle is complementary to the direct measurements of 
CP-asymmetry. It requires mainly studies of the survival probabilities and processes 
where oscillations are averaged or the coherence of the state is lost. 



1 Introduction 



Measurement of CP-violation in leptonic sector is one of the main challenges in particle 
physics, astrophysics and cosmology. 

For three neutrinos (similarly to the quark sector there is a unique complex phase 
in the lepton mixing matrix, 5di which produces observable CP-violating effects 0. (If 
neutrinos are Majorana particles, two additional CP-violating phases exist. These phases, 
the so-called Majorana phases, do not appear in the oscillation patterns.) The phase do 
leads to CP-asymmetry |]^, P(z/a z/^) 7^ P{va — > vp), as well as T-asymmetry P{voi 
Up) 7^ -P(z^/3 i^a), of the oscillation probabilities (see also and references therein). 
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Measurements of the CP- and T- asymmetries provide a direct method of estabhshing 
CP-violation. There are a number of studies of experimental possibilities to measure the 
asymmetries. It was realized that in the 3i/-schemes of neutrino mass and mixing which 
explain the atmospheric and solar neutrino data, the CP-violation and T-violation effects 
are small and it will be difficult to detect them P|. The smallness is due to small values of 
Ue3 (restricted by CHOOZ result) and Am^„„ (responsible for the solar neutrino conversion). 
Still, the effect can be seen in the new generation of the long baseline (LBL) experiments 
provided that the LMA-MSW is the solution of the solar neutrino problem and that 11^3 > 
0.05 

Two types of LBL experiments sensitive to 5d are under consideration |]lOl: the exper- 
iments with superbeams p, Q and neutrino beams from muon storage rings (the neutrino 
factories) Analysis shows ||, |T2| that for {Uesl > 0-05 and Am^^„ = 5 x 10"^ eV^ 

neutrino factories can discriminate between 6o = and 5/) = f at the 3cr level [|T^ while 
according to superbeams are able to distinguish at the 3a level 6^ = from = |. In 
these experiments, the sensitivity to 6d decreases linearly with Am^„„. So, the present un- 
certainty in Aml^^ results in an order of magnitude uncertainty in evaluation of sensitivity 
to 6d in the future neutrino factories and superbeam experiments. If Aml^^ is smaller than 
2 X 10^^ eV^, the direct methods will not be sensitive to 60 [0. Moreover, neutrino factories 
and superbeams are very expensive and technically difficult, interpretation of their results 
can be rather complicated and ambiguous. In view of these difficulties, we need to explore 
any alternative way to search for CP-violation. 

Notice that apart from the asymmetries, the phase So can be determined also from 
measurements of CP-conserving quantities, the oscillation probabilities themselves, which 



depend on 5d [p!4 



The alternative method to establish CP-violation is to measure the area of unitarity 
triangle. This method is well elaborated in the quark sector. Indeed, the area of the unitarity 
triangle, S, is related to the Jarlskog invariant, Jcp, which is a parameterization independent 
measure of CP-violation, as 

S = \jcp. (1) 

So, to establish CP-violation it is sufficient to show that the longest side of the triangle, is 
smaller than the sum of the other two. 

The problem is to measure lengths of the sides of the triangle. As we will see, the 
method of the unitarity triangle differs from measurements of asymmetries and may have 
certain advantages from the experimental point of view. 
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Previously, some general properties of the unitarity triangles for leptonic sector (geometric 
features, test of unitarity) have been discused in [0, [|16], ||T7| . 

In this paper we will consider the possibility to reconstruct the leptonic unitarity triangle. 
In sect. 2, we introduce the leptonic unitarity triangles and study their properties. We 
estimate the accuracy with which the sides of the triangle should be measured to establish 
CP- violation. In sect. 3, we describe a set of oscillation measurements which would in 
principle allow us to reconstruct the triangle. Additional astrophysical measurements which 
would allow us to realize the method are suggested in sect. 4. Discussions and conclusions 
are given in sect. 5. 



2 Leptonic unitarity triangles 

In the three-neutrino schemes the flavor neutrino states, z// = (i/e, t'^, i^r), and the mass 
eigenstates I'mass = {1^1,1^2,^3), are related by the unitary MNS (Maki-Nakag 
0) matrix Q: 







Ue2 


Ue3- 


Umns — 




U^2 








Ur2 


Ur3. 



The unitarity implies 

Ueiu;, + u,2u;, + u,,u;^ = o, 

t/elf/;i + U,2U:^ + f/e3f/;3 = 0, (3) 

Uriu;, + Ur2u;^ + Ur3u;, = 0. 

In the complex plane, each term from the sums in (|]) determines a vector. So, the Eqs. (|^) 
correspond to three unitarity triangles. The CP-violating phase, 6d, vanishes if and only 
if phases of all elements of matrix (Q) are factorizable: Uai = e'^^'^°'^"'^^\Uai\- In this case 
UaiU^i = e'^^'^°'~°''^^\Uai\\Upi\, and therefore the unitarity triangles shrink to segments. 

To construct the unitarity triangle, one needs to measure the absolute values of the 
elements of two rows (or equivalently two columns) in the mixing matrix. The area of 
*The mixing of three flavor states (two light neutrinos and heavy neutral lepton from the third generation) 

have been discussed in jl^ . 
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the triangle is given by the Jarlskog invariant, Jcp Eq. (|I]). The area is non-zero only if 
sin^i) 7^ 0. 



2.1 e — /i triangle; properties 

We will consider the triangle formed by the e- and /i-rows of the matrix @ (see Eq. (^a)). 
(Up to now, there is no direct information about the elements of the third row. Moreover, 
even in future, both creation of intense Vr beams and detection of Vr seem to be difficult.) 
To reconstruct the e — /i triangle three quantities should be determined independently: 

\Ueiu;,i \u,2u;,l \Ue,u;^\. (4) 

The form of the triangle depends on the yet unknown value of |[/e3| and on the specific 
solution of the solar neutrino problem. In what follows, we will consider mainly the LMA- 
MSW solution which provides the best fit for the solar neutrino data. 

In Figs. H and |^, we show examples of the unitarity triangles for different values of Uez 
and 5^1. In these figures we have normalized the sides of the triangles in such a way that the 
length of the first side equals one: 

X = 1, y = ^ and z = ^ . (5) 

\UelU*i\ \UelU*i\ 



We use the standard parameterization of the MNS mixing matrix pO[ in terms of the rotation 
angles 9i2, Ou, O23 and the phase 6d. We take values of 9i2 and 6^23 from the regions allowed 
by the solar and atmospheric neutrino data. 

In Fig. we present the triangles which correspond to sin^ 26'i3 = 0.12 (the upper bound 
from the CHOOZ experiment for Am^^^ = 3 x 10^^ eV^). The arcs show 10% uncertainty 
in measurements of the sides y and z. From Fig. |l|, one can conclude that for maximal 
CP- violation, 5d = 90°, the existence of CP- violation can be established at the 3cr-level or 
even better if the sides of the triangle are measured with 10% accuracy. For 60 = 60°, the 
confidence level is approximately 2a. No statement can be made for < 45° unless the 
accuracy of measurements of the sides will be better. These estimates should be considered 
as tentative ones. In order to make precise statements one needs to perform careful analysis 
taking into account, in particular, correlations of the errors. 

The triangles shrink for smaller values of sin^ 2^13 (Fig. According to Fig. ^ which 
corresponds to sin^ 26'i3 = 0.03, for 6d = 90° CP-violation might be established at ~ 2a 
level. No conclusion can be made for 6d < 70°. 
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The form of the triangle is also sensitive to variations of the angle 612 within the allowed 
LMA region. In Fig. we have set 612 = f and sin^ 26*13 = 0.18. As follows from this 
figure with 10% uncertainty in determination of the sides, CP-violation can be established 
for 5d = 90° and 5d = 60°. 

Note that y ~ 0(1) and z is the smallest side, although its length may not be much 
smaller than others. So, CP-violation implies that 

\Ueiu;,\ < \u,2u;2\ + \Ue3u;,\. (6) 

Similar triangles can be obtained for the LOW and VAC solutions. The unitarity triangle 
is different in the case of the SMA-MSW solution. Taking tan^ 6*12 = tan^ 6sun = 0.0016, 
sin 6*23 = 1/v^ and sin 6*13 = 0.15, we find y = 0.25, z = 0.96. Now y is the smallest side, 
are the two other sides have comparable lengths. Note that in spite of small mixing of the 
electron neutrino the smallest side is not very small. Even in this case a moderate accuracy 
in determination of the sides would allow us to establish CP-violation. 

In general, to establish CP-violation, one needs to construct the triangle without using 
the unitarity conditions. However, if we assume that only three neutrino species take part 
in the mixing and that there are no other sources of CP-violation apart from the MNS- 
matrix, we can use some equalities which follow from unitarity. In particular, we can use 
the independent normalization conditions: 

i=l,2,3 i=l,2,3 

In this case, the CP-violation effect can be mimicked at some level by the 4th (sterile) 
neutrino. To eliminate such a possibility, one should check the normalization conditions 
experiment ally. 

Thus, to find the sides of the triangle we should determine moduli of four mixing matrix 
elements: 

\Ue2l K2I \Ue3l Ksl- (8) 

They immediately determine the second and third sides. The two other elements, \Uei\ and 
|f/^i|, and consequently the first side, can be found from the normalization conditions (|^. 
For the first side we have \U*^Uf,i\ = ^(1 - \Ue2\^ - \Ue3\'^){l - \U^2\^ - |t/^3|2). Taking into 
account this correlation in determination of the sides of the triangle one can estimate accu- 
racy of measurements of the elements (|^) needed to establish CP-violation via the inequality 
(I). Let us introduce 

A = \Ue2\\U,2\ + \Ue3\K3\ - ^J{l - \Ue2V - \Ue3?){l ' P ,2? ' W t^^?) (9) 
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which is a measure of CP violation. CP is conserved if A = 0. For the most optimistic cases, 
where Ue3 is close to the CHOOZ bound and 6d = 90°, we find A = 0.10 - 0.13. 

Suppose the elements \Uai\ are measured with accuracies A|f/aj|. Assuming that the 
errors |Af/cfi| are uncorrelated, we can write the error in the determination of A as 



AA 



, E 1^1 (10) 



where 



dA 


= K2\ 


+ 


\Ue2\\U,l\ 


dA 


d\Ue2\ 


\Uel\ ' 




dA 


= \Ue2 
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\Uel\\U^2\ 


dA 


d\U,2\ 


If^Mll ' 


d\U^3 



\Ue3\ 











u,3\+ ' z (11) 



Ue3\ + ' 'T • (12) 



If^ell 


\Uf.3\ 




U 


Ail 





As an example, let us choose the oscillation parameters used in Fig. 1 and do = 90°. Then 
from Eqs. (0, |l|) we find dA/d\Ue2\ = 0.82, dA/dpe^l = 0-77, dA/d\U^2\ = 2.0 and 



dA/d\U^3\ = 1.9. Note that for muonic elements the derivatives are larger by factor of 2. 
This is a consequence of the appearance of the relatively small element | U^i \ in denominators 
of (0). So, the muonic elements should be measured with the accuracy two times better 
than the electronic elements. 

For our example we find from Eq. ([ToD that AA < 0.065, which would allow to es- 
tablish deviation of A from zero at the 2a level, if A|f/e2| = A|f/e3| < 0.03 and A|[/^2| = 
^If^AtsI < 0.02. This, in turn, requires the following upper bounds for relative accuracies 
of measurements of the matrix elements: 6% for \Ue2\, 17% for |?7e3|, and 3% for \U^2\ and 



\U^3\. Since 



A\U^2\ _ f \U^i\Y A\U^^\ 



\U^2\ \\Ufj,2\ J \U^l\ 

and U^i ~ 0.5f/^2, the required 3% accuracy in \U^2\ corresponds to 12 % uncertainty in U^i. 
If there are correlations between A|f/Qj|, the situation may become better. So, the above 
estimations can be considered as the conservative ones. 



2.2 Present status 

At present, we cannot reconstruct the triangle: knowledge of the mixing matrix is limited 
to the elements of the first row (from the solar neutrino data and CHOOZ/Palo Verde 
experiments) and the third column (from the atmospheric neutrino data). To reconstruct the 
triangle one needs to know at least one element from the block UfSi, where (3 = ix,t, i = 1, 2. 
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That is, one should measure the distribution of the (or/and z/^) in the mass eigenstates 
with spht by the solar Am^. Using the unitarity condition we can estimate only the ranges 
for these matrix elements. Clearly, present data are consistent with any value of the CP- 
violating phase and, in particular, with zero value which corresponds to degenerate triangles. 

Let us summarize our present knowledge of the relevant matrix elements. 

1). The values of the mixing parameters \Uei\ and \Ue2\ can be obtained from studies of 
solar neutrinos. Neglecting small effect due to Ues, for the LMA-MSW solution we obtain 

1^ = I tan I = 0.39 - 0.77, (95%C.L.) (13) 

Well 

and then using the normalization condition: 

\Uei\ ~ [l+tan^ 9 sun]~^^^ = 0.79 - 0.93, (95% C.L.). (14) 



2) . The absolute value of \Ue3\ is restricted from above by the CHOOZ|^ and Palo 
Verde experiments. The 2i/ analysis of the CHOOZ data gives for the best fit value 

lUesl < 0.20, (90% C.L.). (15) 

For lower values of Am^j^, the bound is weaker: |?7e3| < 0.22. 

3) . The admixture of the muon neutrino in the third mass eigenstate, jf/^sl, is determined 
by the atmospheric neutrino data. Again, neglecting effects due to non-zero f/es, we can write 

A\U^3\\1 - \U^3\^) = sin" 29atm, (16) 

where sin^ 29 atm can be extracted, e.g., from analysis of the zenith angle distribution of the 
/i-like events in terms of the — u^- oscillations. Using the Super-Kamiokande data, we find 

\U^3\ = 0.707+°:}2, (90%C.L.). (17) 

4) . At present, there is no direct information about and |f/^2|- To measure these 
elements, one needs to study the oscillations of muon neutrinos driven by Am^^^. The 
normalization condition allows us to impose a bound on a combination of these elements: 

If/^iT + \U^2\^ = 1 - \U^.3\^ = (0-33 - 0.67). (18) 

So, to determine |?7^i| and \U^2\ separately we need to measure a combination of these 
elements which differs from the normalization condition 
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3 Reconstructing the unitarity triangle 



Let us consider the possibility to determine the triangle in the forthcoming and future oscil- 
lation experiments. We suggest a set of oscillation measurements with certain configurations 
(base-lines, neutrino energies and features of detection) which will allow us to measure the 
moduli of the relevant matrix elements (see Eqs. (^, |^)). 

In general, for 3z/-system the oscillation probabilities depend not only on moduli of the 
mixing matrix elements, (^ we are interested in, but also on other mixing parameters includ- 
ing the unknown relative phases of the mixing matrix elements, 5x- Therefore, the problem 
is to select configurations of oscillation measurements for which the dominant effect is deter- 
mined by relevant moduli and corrections which depend on unknown elements and phases 
are negligible or sufficiently small. 

The hierarchy of mass splittings: Am^j„ ^ Am^^„ helps to solve the problem. We use 



e=— ^-0.02 (19) 

as an expansion parameter, where the estimation corresponds to the best fit values of the 
mass squired differences. Another small parameter in the problem is \U(.^\. 

In what follows, we suggest a set of measurements for which the oscillation probabilities 
depend mainly on the relevant moduli: 

PaP = Pap{\Uei\AU^^\) + ^PaP%), a, (3 = 6, ^ , (20) 

where AP ^ P. We estimate corrections, APq,^(5^), due to unknown mixing elements and 
phases. 

It is convenient to study the dynamics of oscillations in the basis of states obtained 
through rotation by the atmospheric mixing angle: (z/g, u'^, z/^), where 

^' = . ^ S ^rSiy, - U,,Vr), K = I ^, A K^^r + U;,V,). (21) 



Projections of these states onto the mass eigenstates equal 



(z.>,) = f/:„ w,w) = - , , Ki^i) = - , (22) 



and 



{ye\y?) = u:,, Kk3) = o, {u'M,) = ^i- We??- (23) 
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Note that in the hmit f/es = 0, the state v'^ coincides with mass eigenstate z/3, whereas 
z/^ = -sin6'i2Z^i + cos6'i2Z/2- 

In matter, the system of three neutrinos {ve-, t'^, ^r) has two resonances associated with 
the two different Am?. The corresponding resonance energies for the typical density in the 
mantle of the Earth, are 

Ef3 = 7GeV, Ef2 = 0.15GeV. (24) 

These energies determine the typical energy scales of the problem as well as the energies of 
possible experiments. Also there are two length scales in the problem which correspond to 
the oscillation lengths: 

These numbers have been obtained for the best fit values of the mass squared differences. 

Let us consider possibilities to determine the moduli of relevant elements of mixing matrix 
(§) in turn. 



3.1 \Ul^U,i\ 

In principle, this product can be directly measured in studies of the z/^ — Vf. oscillations driven 
by Am^^^. Let us consider a relatively short baseline experiment in vacuum. The transition 
probability can be written as 

P,e = m^U,,? sin^ + AP,e, (26) 

where AP^g is the correction due to existence of the Am^^^ splitting : AP^g when 
Amg„„ — > 0. Thus, if the original flux is composed of pure (or pure Ve), detecting the 
appearance of z/g (or z/^), one can measure immediately \U2.^U^-i\ provided that AP^e is small 
enough. Note that AP^e depends on mixing matrix elements Uai-, Ua2, (a = e,/i), both on 
their absolute values and on phases which are unknown. So, we cannot predict AP^g and 
the only way to proceed is to flnd conditions for experiment at which this value is small. An 
alternative method would be independent measurement of If/gsl and |?7^3|. 

For neutrino energies, E > 100 MeV (which are of practical interest) the oscillation length 
in vacuum, I23, is more than several hundred kilometers. This means that the experiment 
should be a long-baseline one, and therefore oscillations will occur in the matter of the Earth. 
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In a medium with constant density^ the probabihty can be written as 



p 

^ fie 



{UTsTu-, (e^*- - 1) + mru;;^, (e^*- _ i) p , (27) 



where f/™ are the mixing matrix elements in matter and $^ is the oscillation phase difference 
of i— and j— eigenstates. 

In the vacuum limit (one may consider a hypothetical configuration of experiment where 
neutrino beam propagates mainly in atmosphere or in a tunnel), = Uai and $™ = 
The first term in (|27|) corresponds to the mode of oscillation we are interested in, and the 
second term is due to the Amf^^ splitting. The main correction follows from the interference 
of these two terms. 

For the correction we find 



U:,U^iU,sU* $32[sin(4 - $32) - sin5j, (28) 



where 6x is the unknown phase of the product of four mixing matrix elements. In derivation 
of (p8|), we have used the smallness of the phase $12: 

$12 = e$32, (29) 

assuming that $32 = 0(1) (which maximizes the effect of oscillations). Then the relative 
correction is of the order of 

AP^e ^ ^ sin26',^„ ^^^^ 

-P/xe I t^eS I 

For the best fit values of the solar oscillation parameters (LMA-MSW solution) and Ues = 0.2 
we get AP^eZ-P/xe ~ 0.1. That is, the product It^es^MSp '^^^ measured with accuracy not 
better than 10% for maximal possible f/gs- Consequently, the accuracy in the determination 
of \U*^Ufj,3\ cannot be better than 5%. 

There are two possibilities to improve the accuracy: 1) The main oscillation term and the 
interference term have different dependences on $32 and therefore on E/L. So, in principle 
one can disentangle these terms by studying the energy dependence of the effect. 2) The 
sign of the interference term can be changed varying E/L. Therefore, the correction can be 
suppressed by averaging over energy, especially if 6x is small. 

Note that for other solutions of the solar neutrino problem (LOW, SMA, VAC), Am^^„ 
is much smaller and the correction is negligible. 



^For simplicity we will consider matter with constant density. Density variation effects do not change our 
conclusions. 
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In the matter the dependence of the oscillation probabilities on mixing matrix elements 
becomes more complicated. However, there are two limits in which the dominant term of 
Pf^e can be reduced approximately to the form (|26|): (i) low energy limit E ^ in which 
matter corrections are small and (ii) short base-line limit L <^l'^^ where "vacuum mimicking" 
condition is satisfied ||23|| . 

Let us consider first the low energy case, E ~ (200 — 500) MeV. The relative corrections 
due to matter effect to the main term in (^) are of the order of 

/23 2^/2GFneE 



Iq Am, 



(31) 

atm 



where Iq = \/2'n /Gpne is the refraction length. For E ~ 200 MeV, we have e ~ 0.02, while 
for ~ 1 GeV, the corrections reach 10%. Moreover, the matter effect is of order 1 for the 
correction term driven by Am^^„. 

At low energies the mixing in the heaviest eigenstate is only weakly affected by matter, 
so that in the first approximation we can take ~ Uez-, UJ^ ~ f/^s. For E ~ 200 MeV the 
oscillation length due to 2 — 3 level splitting is ~ 200 km, and therefore the optimal baseline 
would be L ~ (100 - 200) km. 

The energies E ~ 200 MeV are in the resonance interval for Am^^„. This means that the 
electron neutrino has comparable admixtures in the two light eigenstates: ~ ~ 1/ ^/2. 
The oscillation length is of the order of the vacuum oscillation length (for the LMA-MSW 
solution /i2 ~ 5 ■ 10^ km), and therefore the oscillation phase due to the 1-2 level splitting is 
small: $^ ~ ItiL / 1\2 ~ f ^ 1- These features simplify the analysis of the correction term. 
Indeed, using (|2l|, |2|, |2|)), we find 



Therefore, the relative correction AP^eZ-f^te which appears due to the interference of the 
term (^) with {U*^)'^UJ]^ in (^ can be written as 

AP^e sm29sun\Ur3\ slu 26'^„„ 

~ ^ — i i — I \ ^ ^ — i i — ■ ^ ) 

I ^e3 1 I I I ^e3 1 

Here we have taken into account the relation between phases (^) . So, the expression for 
the correction is basically reduced to that in the vacuum oscillation case given by Eq. (|30D. 

Similar considerations hold for the antineutrino channel. We find that not only the main 
term but also the orders of magnitude of the corrections coincide with those for the vacuum 
case. 
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Let us consider the case of high energies, E ~ (5 — 10) GeV, and relatively short baselines, 
L ~ 700 km, for which the "vacuum mimicking" condition, L <^ /J^, is satisfied. Notice that 
these energies are in the range of the resonance due to Aml^^ (p4|). 

In the limit where the Am^^„ splitting can be neglected the probability (for the constant 
density case) becomes 

P,e = Sin^ l$32i?i3 + AP,e, (34) 



where 

RL = I cns2/9io - 



i??3 = (cos 2^13 - ^ W sin^ 2^13. (35) 



The oscillation phase in matter, can be written in terms of the phase in vacuum, $32, 
as 

$^2 = '^'32i?13, $32 = 2nL/h2. (36) 



For <l>^ <C 1 the expansion of the probability (|54D in powers of $^ leads to 

P,e = MU*M'<^lJl-^<!>lX (37) 



Note that the first (leading) term in ( |37D reproduces the vacuum probability (vacuum mim- 
icking) . 

The correction depends on energy. In the resonance R^^ = sin^ 26'i3, and therefore accord- 
ing to (|37|) corrections are below 10%, even if $32 ~ 1. Below the resonance: for /23//0 0, 
we have R^^ 1 and to have a small correction $32 should be small, thus leading to sup- 
pression of the oscillation effect. Therefore, for such a type of experiment, the optimal range 
of energies is = (5 — 10) GeV, and the optimal baseline ~ 700 km. Above the resonance 
where Rf^ 00, the oscillation effect is even more suppressed. 

Let us evaluate the correction to (^), AP^e, due to the solar mass splitting. For high 
energies both oscillation phases $^ and $^ are small, so that from Eq. (^) we obtain the 
interference term 

AP,e = 2Re [(t/S)*f/,l$3"^2f/e^(t^ri)*'^r2] • (38) 
Considering explicitly the mixing in the — v'^ system we find 



myu^3^32 ^ u:,u^,^,,. m 

The mixing of the electron neutrino in the first eigenstate is strongly suppressed by matter 
effect: 

[/-([/-)* ^ U,,wJ^ , (40) 

tl2 
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where I0/I12 ^ 1- As follows from the level crossing scheme, in the resonance region the 
phase difference between the two light eigenstates is 



~ -$32. (41) 

Plugging the expressions (0), (0) and (^l]) into Eq. (|38D we find 



U:,U,,U,,U;,] ^l,^. (42) 

tl2 



The relative correction can then be written as 

AP^^ sin2^^ 

Near the resonance the relative correction is similar to that in the low energy limit or in 
vacuum case. The correction can be further suppressed if E' > E^^. At the same time, with 
the increase of energy, the phase of oscillations ($32 in Eq. (pTj)) decreases, and therefore 
the number of events decreases quadratically. 

Let us consider the antineutrino channel. At high energies, ~ ^'^ and z/g mixes with 
u'^ in the uim and eigenstates. In the limit i/2m = t'^j, the standard vacuum expression 
for the probability Pe^ (Eq. (^) with R13 R13) is reproduced. 

The correction is related to the admixture of z/g in U2m which is determined by sin 6^ and 
is strongly suppressed by matter. A straightforward calculation give 

Ki3 K12 COS ^13 

where the correction is given by the second term. Taking into account that $'5^ = $i2-Ri2 = 
— $32-Ri2e one can find 

sin26'5„„ 





sin 29sun 


|f^r3| 


COS^^^ 


P 

fie 


~ e 

f/e3 


\Un3 


COS 6'i3 



Note that in contrast to the neutrino case the correction does not change with energy. So, 
using neutrino beam seems to be more promising because for neutrinos relative corrections 
decrease with increase of energy. 

We now discuss the sensitivity of upcoming and planned experiments to the product 
l^es^AtsI- shown that combining the data from the MINOS and ICARUS ex- 



periments, one can obtain an upper bound sin^ 26*13 < 0.01 at the 95% C. L. This would 
correspond to [U^.^U^'^l < (0.03 — 0.04) which is about 4 times stronger than the present 
bound: |?/e3^M3| < 0.15. The searches for the z/^ — Uf. oscillation will be performed in phase I 
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of the JHF project. The sensitivity to the product \U*^Ufj_:i\ can reach to 0.02 Therefore, 
if \U*^Uf^3\ is at the border of the present upper bound it will be measured with about 15 
% accuracy. Neutrino factories will be sensitive to It^es^Msl down to 10~^. However, for 
\U*^U^3\ < few X 10~^ the correction due to non-zero value of Am^„„ (in the case of LMA- 



MSW solution) will be comparable with the main term (see for related discussion). For 
other solutions the corrections are negligible. 



3.2 \Ue3\ 

Independent determination of {Uesl seems to be important in view of the difficulties associated 
with the direct measurements of \U*^Uf^3\ discussed in sect. 3.1. Knowledge of If/gsl is also 
needed for a precise determination of \Uei\, |t/e2| and other mixing elements. 
The survival probability for i/e-oscillations in vacuum can be written as 

2 



|f/e3p (e^*- - l) + 1 + |t/eiP (e^*" - l)| . (46) 



Note that, in contrast to the conversion case, the probability amplitude depends on the 
required moduli of the matrix elements. A similar analysis holds for antineutrinos. 

For low (reactor) energy experiments the matter effects are negligible and the probability 
equals 

Pee = 1 - 4(1 - \U,s\')\Ue3\' sin^ ^ + APee . (47) 
Here the correction APee due to the Am^„„ splitting can be evaluated as 

AP,e = 2|t/,in[/,3|2$,2sin$32 " slu^ 2^,,„. (48) 

The relative correction is small: APee/ Pee < 2 %, so that in principle, jf/esl can be determined 
with better than 1% accuracy. Experimental errors in the measurement of Pee will dominate. 

Let us comment on the experimental prospects for measuring jf/esl- A new reactor ex- 
periment, Kr2Det, has been proposed which will be able to set the bound jf/esl < 0.07 at the 



90 % C. L. ||2^. This bound can be used to estimate the sensitivity. If, e.g., \Ue3\ = 0.2, one 
would expect that the experiment will give |f/e3p = 0.040 ±0.005. Consequently, jf/esl itself 
will be determined with about 6% accuracy. In fact, the situation can be slightly better. If 
lUesl is near its upper bound, one can study spectrum distortion and therefore to perform a 
more accurate determination of {Uesl- 

It is not clear if future measurements allow us to measure \Ue3\ precisely enough to 
reconstruct the third side of the triangle. But certainly, they will contribute to a more 
precise determination of \Uei\ and \Ue2\- 
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3.3 lU^sl 

The present analysis of the atmospheric neutrino data in terms of 2z/- mixing gives the bound 
(p!7|) on If/^sl- Note that at this accuracy, the effects of non-zero Ues and subleading modes 
driven by Am^„^ are unimportant. Indeed, in P3, it is shown that the allowed region for 
If/^sl does not change considerably as {Uesl varies between zero and its maximal possible 
value. To reconstruct the unitarity triangle, we need a more precise measurement for |t/^3|, 
which requires 3z/ analysis taking into account the effect of non-zero |?7e3|- 

The element |f/^3| can be measured in z/^-disappearance due to oscillations driven by 
Am^^^. The z/^-survival probability in a uniform medium equals: 



P 

J- tifj. 



prn^2 l^gi*- - l) + 1 + |f/-|2 (g'^r^ _ l) |' . (49) 



Again, there are two limits in which the dominant term of this probability reduces to the 
vacuum oscillation probability plus small corrections: (i) the low energy limit {E <^ -E^), (ii) 
and the high energy case ( E > Er) with a small baseline for which the vacuum mimicking 
condition is satisfied. 

Let us consider first the high energy limit. The dynamics is particularly simple, ii E > 
E^^. In this case ~ t'^, whereas the states and u'^ strongly mix in z/2m and u^rn- For 
the phases we have ~ $32 and ~ $32-^13 ^ 1- Neglecting the admixture of z/^ in 
uim and i'2m (which is smaller than |f/e3|e) we obtain 



\TT 12 \Tr 12 
TTm\2 I^r3| |r,m|2 I^m3| 2 /c-qN 



where cos^ 9^ = [1 + (cos2^^i3 — Ii^/Iq)/ Ri3]/2. Inserting these matrix elements into (^) we 
can reduce the probability to the form 

P^^ = 1 - (1 - \U,s\')\U,s\'e,, + AP,, , (51) 

where AP^^ is the correction due to matter effects, non-zero value of Ue3 and Am^„„. Re- 
calling $^ <^ 1, we can write 



f/e3rif/^3r -2|f/^3r(i - lu^sniRncos'e^^, - lu^sn n^- (52) 



At the resonance, P13 cos^^ ~ Wes] and the corrections are strongly suppressed. However 
the resonance region is rather narrow. Above the resonance, P13 cos^ 6*]^ ^ |f/e3p/Pi3, and 
as follows from (|52|), the relative corrections can be estimated as AP^^/P^^ ~ |f^e3p < 0.04. 
Moreover, the corrections are calculable and can be taken into account once |f/e3p is measured 
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even with a reasonable accuracy. Corrections, which depend on f/^i and ?7^2, are suppressed 
by the ratio e at ~ E^^. 

In the antineutrino channel we find 



\U. 



m |2 





U,3 


2 


1 - 




u 


e3 


2 



1 - 



eSl 



ai3 



(53) 



where ais ^ 1 - 2|f/e3p + /13//0, '^^2 ~ "^32, and 



Irrm 1 2 
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f/e3 
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f/r3 


sin 26'5„„ 




1 - 
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e3 
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«13 




-Rl2 


ai3 



COsSr 



(54) 



Here (5a: is the unknown phase of the product of matrix elements. We can estimate $ 



m 
12 



^13 



^23 ~ 



$^ + $32 ~ $32(1 - -R13). Using the first term in the Eq. 



lowest order terms in |f/e3p, we find 



MM 



1 - «i2 



((1 



m3 I 



|t/^3r + |f/e3rif/^3r(i-ar3)(i + « 



'^13 



2\a 



and keeping the 
M3l))- (55) 



As in the neutrino case, we assumed $^ <^ 1. According to (0) the correction to the 
standard expression for 2z/ probability (the first term in the brackets) is of the order of |f/e3p 
with coefficient smaller than 1. Note that, the unknown phases are not involved in the Eq. 
(^), so the corrections will be calculable once {Uesl is measured. 

The relative corrections which depend on unknown phases originate from the interference 
of the term proportional to |f/^P (the third term in Eq. (^91) replacing $™ 
\UJ2\) and the main term. They can be estimated as 



$^ and 



\^ fill 



AP, 



fj.fl 



I -P. 



MM 



2esin26'5„„|f/e3|. (56) 

Note that in contrast to the case of determination of \U*^U^'^\, the relative corrections are 
suppressed by \Ue3\ because, in this case, the main term is larger; it corresponds to the 
dominant mode of oscillations {i.e., P^^ ~ 1, while P^e ^ !)• 

Let us consider the low energy experiment with E ~ E^2 ~ (200 — 500) MeV. In this case 
the Am^j^-driven oscillations are in the quasi-vacuum regime {UJ^ ~ f/^3, $^ ~ $32) and 
the base-line can be relatively small: L ~ /12 ~ 100 km. On the other hand, the oscillations 



driven by Am^^^ are in the vacuum mimicking regime: $^ ^ 1. It can be shown that, 



I rr™ 1 2 



\Ur3 



-sin^^^-. 



^12 ~ 



i?12$l2 = e-Rl2*^'l3, 



(57) 



where sin^ 9^ 
2) system: 



1- I^e3, 

[1 — (cos 26*12 — ^i2/^o)/-Ri2]/2, and Ru is the resonance factor for the (1 



R 



12 



cos 2^12 -y^l +sin^2^i2. 



(58) 
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Inserting the matrix element ( pTD into (^), we can reduce the probabihty to the form of Eq. 
(IH) with 

AP^^ = 2e|f/.3r|t/,.3|'$32sin$32Pi2sin2^r2 • (59) 

Let us consider last two factors in this expression. In the resonance, Pi2 sin^ ^5"2 ^ P12/2 = 
sin2^s^j„/2, but above it Ri2sm^9^ lu/^o and the correction increases with energy too. 
Below the resonance R12 sin^ sin ^o„„. 

Thus, in the resonance region and below it, the correction is small and of the order of 
^"^Ln/^"^atm- The corrcctious due to admixture of u'^ in the lowest mass eigenstate (which 
we have neglected) are of the order |f/e3p. As in the high energy limit, the relative corrections 
are restricted to AP^^/Pf^^ < 0.04, and moreover, the dominant part of these corrections 
can be calculated in terms of jf/esl- 

In the antineutrino channel, similar consideration gives the following corrections which 
can be calculated in terms of the moduli of the matrix elements: 



|2 



^^M/. = ^ !r'\!r ^2 ^32 Sin $32(i?12 ' \l ^2 ' sin' 2esun). (60) 

The relative corrections are of the order of Am^„„/Am^j„. The corrections which depend 
on unknown phases, are further suppressed (~ |f/e3|e)- 

Studying the disappearance of z/^, the MINOS experiment will determine Am^^„ and 
(1 — |?7^3p)|?7^3p with 10 % accuracy at the 99% CL. after 10 kton-years of data taking 
p9|, |2^. Much higher precision can be achieved in phase I of JHF: the oscillation parameters 
((1 — |?7^3p)|?7^3p and Am^^,^) will be determined with 1% uncertainty 0. Thus, there are 
good perspectives to determine |f/^3| with precision better than 2 - 4 %. 

Notice that the future atmospheric neutrino experiment, MONOLITH, can measure 
sin^ 26^23 with uncertainty of 8% ||30l . 



3.4 I [/el I and \Ue2\ 

The values of \Uei\ and \Ue2\ can be obtained from the solar neutrino data. To first approx- 
imation, due to the low energies of solar neutrinos the matter effect on \U(.z\ is negligible 
and the solar neutrino conversion driven by /S.m?^^,^ will produce only an averaged oscillation 
effect. In this case the survival probability equals |^ 

Pee = (1 - |?7e3r)'P2(tan2 Osun. Am^^J + \U,^\\ (61) 
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where 



tan^ Osun = TTT^ (62) 

Well 

and P2 is the two neutrino oscillation (survival) probability determined from the solution 
of the two neutrino (z/g — u^) evolution equation with the oscillation parameters tan^ Ogun, 
Aml^^ and the effective potential (1 — If/esHK- 

Precise measurements of |t/e3p will be performed by the KamLAND experiment for which 
Pee is given by ( |6TD with P2 determined by the oscillation formula in vacuum: 

"2 = 1 — 7 ; r^TTTT sm . 63 

The expected error in determination of sm'^29sun, and therefore the combination |t/eiP|f^e2p 
is around 5% Then using the measured value of |?7ei||f^e2| and the normalization con- 
dition, |?7eiP + |f^e2p = 1 — l^esP, wc cau find |?7ei| and |f/e2|, separately. The accuracy can 
be better than (2 - 3)% . 



3.5 \U^i\ and \Uf,2\ 

The determination of |f/^i| and |f//^2| is the most challenging part of the method. Note that in 
contrast to \U*^Ue^\ (see sect. 3.1), it is not possible to measure the combinations \U*-^^U^i\ or 
\U*2U^2\, directly from the oscillation experiments. Indeed, in vacuum the z/^ — Ue transition 
probability is determined by the product Re U*iUeiU^2U*2 which depends not only on 
the absolute values of the matrix elements but also on their phases. (For example, in the 
case that Am'j.^^L/E is not resolved, the probability Pe^ is determined by the combination 
W^iUei + U*2Ue2\-) Therefore we will consider the possibility to measure separately and 
\U^2\, so that the second side of the triangle can be found using the electron matrix element 
\Ue2\ obtained in other experiments. In fact, it is sufficient to measure some combination 
of \Uf^i\ and \U^2\ which differs from the normalization condition (p^S)). This requires an 
experiment sensitive to the splitting between the first and second levels associated with 
Amj.^^ which appears usually as a subdominant mode. To suppress the leading effect and the 
interference of the leading and sub-leading modes, the oscillations driven by Am^j^ should 
be averaged out. This condition necessitates the following experimental configuration: 

1) . The energy of beam should be low: E < 1 GeV. 

2) . The baseline should be large: L ^ I23 (in contrast to configurations considered in the 
previous subsections). Moreover, to avoid suppression of the subdominant mode we need L 
to be of the order of the oscillation length due to the (1-2) splitting. 
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At E < 0.5 GeV, we have I23 ~ 500 km, and consequently, to reach averaging the basehne 
can be ~ 2000 km. In this case ~ 0(1). 

To produce muons, we need E > 100 GeV. For these energies matter effects on (1 - 2) 
mixing are non-neghgible and moreover, since the basehne is large, no vacuum mimicking 
will occur. 

The experiment we have arrived at, seems even more difficult than that for direct mea- 
surements of the CP-asymmetries p|. However, our proposed experiment measures quanti- 
ties different from asymmetries, and moreover, only one beam, neutrino or antineutrino, is 
sufficient. 

Let us consider the — oscillation [disappearance] experiment with E ~ E^2 ~ 
(200 — 500) MeV and L > 2000 km. At these energies the inffuence of the matter effect on 
flavor mixing in the third mass eigenstate is small so that we can take ~ Ues and also 
UJ^ ~ Ufi3- (The corrections are of order of e.) Therefore, the normalization condition gives 
+ I^SP — ^ ~ l^esP ~ 1 — |?7e3p- The mixing is reduced to the mixing in 2i/-system, 
so that matrix elements in matter can be obtained by substituting Uei — ^ U^, Ue2 U^. 

The general form of the probability in a medium with constant density is given by Eq. 
(pl|). Let us calculate \UJJ^\. First we express the vacuum value of U^i in terms of Uei, 11^2 
and mixing in the third state, (a = e, /i). To do this, we use U^i = Sq,(z/^|z/^)(z/^|i^i) 
and the relations (^, ^ |2^). Then in the expression for U^i = U^i{\Uei\, \ Ue2\, lUasl) we 
substitute Uei ~^ U'^ (z = 1, 2). A straightforward calculation gives 

= Ir ,2 [l(^e")*t^.3| + U2U:,U,:^ . (64) 

^ " |l^e3| 

Mixing elements in matter can be written as 

lrrm|2_|rr 12 ^03^^!^ _ |^el P ^12 + COS 2^12 - /12//0 , . 

I^eil -l^^il eos^^is " R12 l + cos2^i2 ^ ^ 

(here cos 2^12 = 2|?7ei|V(l - If/eaH - 1), 

\Ur2\' = ^ - lUTll' - \Ue3\' (66) 

and 

{urirur2 ^ ^u:,Ue2. m 
-n.12 

Using these equations we can express in terms of the mixing parameters in vacuum as 

\Ki\' = ^\U,i? + F, (68) 
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where 

^ = T^ \\Ur,?U + |f/e3r|f/M3|V-l (69) 

U — |<^e3| ) '- 

and 

i?i2 - 1 ± /12//0 

- — m;, — • 

Note that in the vacuum hmit /± ^ 0, R12 — > 1 and F — >■ 0. At the resonance, R12 sin 29 12 
and above the resonance where E ^ i?^^ 

iT-r 12 

In this case the dependence of |f/^P, and consequently of the probabihty, on \U^i \ disappears 
in agreement with our result for the high energy version of the experiment in sect. 3.3. The 
survival probability can be written as 

P.. - if/.3r + (1 - \u,.?f - 4 + ^) (1 - i^-^i' - ^ - ^) ^ • (^2) 

Let us underline that F = F{\Uei\'^, l^^asP) is a known function of |?7eiP and If/asp and it can 
be determined once these elements are measured. The contribution of the |?7^i [-dependent 
terms to the probability is about 10%. Therefore to determine precisely enough, the 

probability should be measured with better that 1% accuracy. 

The correction to the formula ( |72D due to matter effects are of the order e. 

For antineutrinos the probability is given by expression (|72D substituting, I12/I0 
— /12//0, R12 — ^ -R12, ^r2 ~^ (obviously, \Uai\ = \Uai\)- Note that in this case, above 
the resonance {E > E12) we get 



and again the dependence on \U^i\ disappears. 

In general the aforementioned conditions (to measure |f/^i|) are fulfilled for the sub-GeV 
atmospheric neutrinos reaching the detector through nadir angles between 30° (for which 
the baseline is tangent to the core) and 80° (with L ~ 2000 km). Indeed, for such neutrinos 
the phase of oscillations driven by Am^^^ is of order 1: Am'j.^^L/2E ~ VeL ~ 0(1), while 
Aml^^L/2E ^ 1. However, due to the presence of both electron and muon neutrinos in the 
initial flux, the number of observable events, e. g. /i-like events, depends both on survival 
and on the conversion probabilities {Pefj. and -P^e)- One can easily show that for conversion 
probabilities, the effects of interference terms, which depend on unknown phases, are non- 
negligible. So, it is not clear, whether atmospheric neutrino data can help to measure \Uf^2\- 
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4 Do alternative methods exist? 



A straightforward (and similar to what we do in quark sector) way to determine the ele- 
ments of the MNS matrix (and therefore the sides of the unitarity triangle) is to study the 
charged current interactions of neutrino mass eigenstates, ui. Indeed, the cross-section of 
the interaction 

i/, + X ^ Z + y, 

where I is a charged lepton, is proportional to |C/hP. In particular, measuring the number of 
electrons and muons produced by the z/i-bcam one can immediately find the ratio |C/ei|/|t^Mi|- 
To perform such a measurement one needs to create a beam of pure neutrino mass eigenstate 
energetic enough to produce the charged lepton, /. There are several ways to produce (in 
principle) a pure mass eigenstate beam: (i) via adiabatic conversion, (ii) due to spread of the 
wave packets and (iii) as a consequence of neutrino decay. In general, one can also use a beam 
of several mass eigenstates provided that they are incoherent. Processes induced by such a 
beam will be determined by the moduli of matrix elements. Effective loss of coherence occurs 
due to averaging of oscillation of neutrinos from far objects (for which Aml^^L/2E ^ 1). 
We will consider these possibilities in turn. 

4.1 Adiabatic conversion of neutrinos in matter 

In a medium with high density (larger or much larger than the resonance density) mixing 
can be suppressed. That is, the flavor state, produced at such a density, coincides with the 
eigenstate of the instantaneous Hamiltonian: Uf ^ z/j^- If the density decreases slowly to zero 
along the path of neutrino, such that the adiabaticity condition is fulflUed, the neutrino state 
will always coincide with the same eigenstate: z/(t) ~ i^imit)- As a result, when the neutrino 
exits the layer (at zero density), it will coincide with the mass eigenstate i'{tf ) ^ = Vi- 

This happens for solar neutrinos with energies 5-14 MeV in the case of the LMA-MSW 
solution. The electron neutrinos produced in the center of the Sun are converted to z/2-state. 
So, by studying the interactions of neutrinos from the Sun we can measure |f/e2|- 

Obviously, usual solar neutrinos cannot produce muons. Measurements of \U^i\ and/or 
\U^2\ will be possible, if high energy neutrinos {E > m^) appear in the center of the Sun 
and propagate adiabatically to the surface. Such a possibility can be realized if massive dark 
matter particles, WlMPs, are trapped inside the Sun and annihilate emitting neutrinos. 

Suppose that the dark matter is composed of neutralinos, The neutralinos annihilate 
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into the Standard Model particles: XX ~^ W+W^, ZZ, qq etc., which in turn decay producing 
neutrinos and antineutrinos. The energy spectra and the flavor composition of neutrino 
fluxes (as well as the absolute value of the flux) depend on the parameters of the SUSY 
model. Generically, one expects an asymmetric flavor composition. Indeed, neutralinos 
annihilate preferably into 66, rf (and if they are massive enough also into tt, W^W~ and 
ZZ). Moreover, muons, pions and kaons are absorbed or lose a substantial fraction of their 
energy before decay. In contrast, the r-leptons decay before appreciable energy loss. As a 
result, one expects an excess of Vr and approximately equal fluxes of and v^: 

F\ur) = F\ur) > F\u,) = F\u,) ^ F\u,) = F\u,). (74) 

At high energies (-E > 100 GeV) the inelastic interactions inside the Sun are very important 
and due to differences in the cross-sections one expects different energy spectra for neutrinos 
and antineutrinos. However, for E^, < 50 GeV the effect of inelastic interaction is smaller 
than 10%. 

Note that in contrast to the case of usual solar neutrinos (for which pure z/g flux is 
generated), WIMPs produce a neutrino flux with a complex fiavoi composition. This creates 
two problems: (i) one needs to know the flavor composition which is subject to various 
uncertainties; (ii) the flnal flux is a mixture of mass eigenstates (and is not a pure mass 
eigenstate). 

Another problem is that only for rather low energies, the adiabaticity condition is fulflUed 
in the resonance channel. For E > 1 GeV neutrinos cross two resonance regions inside the 
sun: the high density (/i)-resonance associated with Am^^^ at density p < 30 g/cm^ and 
the low density (/)-resonance associated with Am^„„ at p < 0.5g/cm^ for the LMA-MSW 
solution. For deflniteness we will consider the scheme with normal mass hierarchy in which 
both resonances are in the neutrino channels. 

The jump probability at the resonance which characterizes the adiabaticity violation can 
be written as Pc ~ exp(— 7 sin^ 9), where 

/ Am2 \ /lGeV\ , , 



where we have used the density proflle of the Sun in ||3^ . (The above formula is valid only for 
weak violation of adiabaticity: ^ 1.) For E = A GeV in the 1-resonance associated with 
Am^„„ = 5 • 10~^ eV^ and tan^ 6 = 0.35 we obtain Pc ~ 0.01. At the h-resonance violation 
of adiabaticity is negligible, provided that Ugs is not very small. We have Pc ~ 10^^, for 
It/esP = 0.03. 
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So, the adiabaticity violation effects are below 5 % for E < 5 GeV. In the antineutrino 
(non-resonant) channel the adiabaticity violation occurs at larger energies. 

The effects of adiabaticity violation lead to the appearance of interference terms which 
depend on unknown complex phases. Therefore one needs to select low energy events. On 
the other hand, even for light neutralinos, ~ 50 GeV, we expect that only about 30% 
(or less) of neutrinos have energies less than 5 GeV. This diminishes statistics and hence the 
efficiency of the method. Notice also that large underwater (ice) detectors have rather high 
energy thresholds, so detection of GeV neutrinos is may be problematic. 

In what follows we present a simplified consideration (neglecting the inelastic interactions) 
to illustrate possibilities of the method and its shortcomings. A detailed analysis will be given 
elsewhere p5[] . 

Using relations (0) we can write the 1^/3 flux [jS = e, fi, r) at the Earth as 

= E = + ^F'Ppr, (76) 

a 

where F° = F° = Fj^ is the common flux of the electron and muon neutrinos, AF" = F^ — F^ 
and Pj3a is the z/q, Z//3 conversion probability on the way from the production region in the 
center of the Sun to the detector on the Earth. (Here for simplicity we do not consider the 
Earth matter effect.) A similar expression can be written for the antineutrino channels. 

Neutrinos from WIMP annihilation can be detected by large underwater and ice detec- 
tors via charged current interactions. In these detectors the rates of /x-like events will be 
measured. The detectors will not be able to identify the charge of the produced lepton and 
therefore we need to sum the neutrino and antineutrino fluxes in our analysis. Using Eq. 
([76|) we can write the expression for the rate of the /x-like events as 

= K + J dn^F\Pr^a^ + P,^a^), (77) 

where 

Nl = jdnF\a, + a,) (78) 

is the rate without oscillations. In the above equations, o"^ and are the cross-sections 
of the charged current interactions of neutrinos and antineutrinos, respectively. Here / dVt 
includes the integration over the neutrino energy, the angle between the neutrino and the 
produced muon and the energy of muon. One should also include the efficiency of detection 
and the energy resolution function. 

Let us find the transition probability, P^r-, which determines according to ([76|) the oscil- 
lation effects. The general expressions for the probabilities Ppa are given in Ref. p6[. Here 
for illustrative purposes we will consider the case of pure adiabatic propagation in the Sun. 
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For E < 5 GeV and the relevant Am^ the oscillatory terms will be averaged out, in 
particular, due to finite energy resolution of the detector and change of the distance between 
the Sun and the Earth. Taking into account this averaging effect we find the Ur vp 
conversion probability in the adiabatic limit: 

Ppr= E \Up^?\U':^^\ (79) 
i=l,2,3 

where f/™ is the mixing parameter in matter at the production region. 

The density in the production region is much higher than the resonance densities, so that 
the mixing is strongly suppressed. Considering the level crossing diagram, it is easy to show 
that u^rn ~ —Ve, i^2m ~ K ^im ~ — t'^- From thcse relations and the definition of z/^ and 
K dHI) *6 obtain 

\U ql ir/ ql 

Inserting (^Oj) into (|79D we find the probability of the conversion: 

1 

^e3| 

Since the atmospheric mixing is close to maximal: |f/^3p ~ 1/2, the dependence of the 
probability Prfj, on |?7^2p is weak. 

In the antineutrino channel, at high densities we have vim ~ z^e, ^2m ~ ^'^ and v^m ~ ^'r- 
Consequently, the mixing elements equal: 

1/7™ I ^ n 1/7™ I ^ l^^^l 1/7™ I ~ 



p = . 

1 - It/. 



U,:>,\\l-\U^-i\') + {l-2\U,^\')\U, 



m2| 



^11 



Using (|7^) we find the probability of the u^. — > z/^ conversion: 

= -r^^nrn V " l^/^sl' - \Uez? + ■ (83) 

Here |f//^2p appears with a relatively large coefficient. 

In the adiabatic limit, the conversion probabilities do not depend on energy and the 
expressions for the rate of events can be written as 
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U,3 




1 - 


Ue3 
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^i2\ 



J dQAF'lil - 2\U,,\^)a^ + {U.^l^cr^]. (84) 
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According to (|5^, the relative effect of the term proportional to |f/^2| is suppressed by 
smaller value of the antineutrino cross-section (t^/ct^ ~ 1/2. 

The relative contribution to number of events from the term which depends on |f/^2p at 
low energies (see (|8^)) is 

^ ^ '^"'^0^"''' / dnAF'a,. (85) 

For larger energies E > 7 GeV, in the neutrino channel the effects of the adiabaticity violation 
are ~ 10%, i. e. larger than the level of required accuracy in the determination of the mixing 
elements. In the antineutrino channel the adiabaticity violation is weaker. So if detector 
is able to identify the charge of lepton, and consequently, to select antineutrino events, one 
will be able to perform better measurements. In particular, events with higher energies can 
be studied. 

Let us comment on the possibility to detect neutrinos from WIMP annihilation and to 
measure |f/;i2p- The event rates due to these neutrinos in a km'^-size detector can be as 
large as few 10^ events/year [^], although the rate is very model dependent. If AF^/F^ ~ 
0.2 — 0.5, the contribution of the term sensitive to |f/^2p is about 10 %. Therefore |t//i2p 
can be determined with accuracy 10% at best, provided that all other involved parameters 
are known. In particular, one should know the original flux and the difference of fluxes 
F° — F° as functions of energy. 

There are several possible ways to deduce information about the ratio of original fluxes 

1) . Theoretical predictions: In principle, future high energy experiments at colliders {e.g. 
LHC) , as well as results of the direct searches for dark matter particles will help to measure 
the mass and the composition (Higgsino-like versus gaugino-like) of neutralinos. This will 
allow to predict the relative neutrino fluxes from annihilation. 

2) . Information on relative neutrino fluxes from WIMP annihilation can be obtained by 
detecting neutrinos from WIMP annihilation in the Earth center. 

These studies cannot determine the absolute value of the original flux [F^). Once we 
obtained by the aforementioned methods the value of AF^/F^, we can try to measure both 
|?7^ip and the original fluxes from studies of solar neutrinos themselves. If the detector is 



able to identify the flavor |^8[|, we can compare the rates of /x-like with r-like or e-like events 
to find the total flux and the value of If/^^ip. 
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4.2 Spread of the wave packets 

Bunches of neutrino mass eigenstates can be obtained as a resuh of difference in the group 
velocities. Neutrinos with a mass squared difference Am^ but the same energy E, produced 
in a source at the same time, will arrive at the detector with a time difference 

/ L \ f Am^ \ /lOO MeV\2 , , 

At = 0.1sec — ^ . 86 

V1028 cm/^SxlO-s eVV V E J ^ ' 

Here L is the distance from the source. If the time during which neutrinos are produced 

at the source, Tp, is considerably smaller than At {vp < At), and if the energy spread is 

small enough (or the detector is able to select neutrinos of certain energy), neutrinos will 

arrive at the detector in bunches: the heavier neutrinos arrive after the light ones. We 

can measure the numbers of charged leptons produced by different bunches via the charged 

current interactions. Thus, the ratio of number of muons and r-leptons produced by the 

first bunch gives \Uni/Uri\. Similarly, the second and third bunches give information about 

\Ui2\ and \Ui3\, respectively. The number of charged leptons of a given flavor, /, produced by 

the first and second bunches is proportional to |f/zi|/|f/z2|5 etc. 

According to (0), the time difference in arrival of the bunches for Am^ = 3 x 10^'^ eV^, 

E = 100 MeV and L = 10"^^ cm equals At = 0.1 sec. So, the duration of the neutrino pulse 

should be smaller than 0.1 second. Moreover, the number of events induced by a single pulse 

should be large enough. It is not clear if the required sources of neutrinos exist. 

4.3 The decay of neutrinos 

The neutrino decay provides another possibility to get pure beam of mass eigenstates. In 
the minimal extension of the Standard Model (in which neutrinos are massive and there 
are right-handed neutrinos) neutrinos can decay radiatively: uj ^ z/j + 7. However, the life 
time is extremely large: r^, > 10^^ s. In certain extensions of the SM the radiative decay or 
3z/-decay may be much faster, however, according to astrophysical bounds, the lifetime of 
radiative decay must be much larger than the age of the Universe (see review [^]). 

The decay which satisfies all the bounds and is relevant for our analysis is the Majoron 
decay 11,11: 

Ui Uj + J, (87) 

where Ui and uj are mass eigenstates, and J is the Majoron. 

Let us assume that t^, ^ 10^'^ sec, so the neutrinos from the Sun do not decay and the 
solar and atmospheric anomalies are explained by oscillations while neutrinos from very far 
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sources {i.e., the gamma-ray bursters, the Active Galactic Nuclei and supernovae) decay 
before reaching the Earth. Then at the detectors the neutrino flux from the far source is 
composed only of the lightest neutrinos ui and Ui. 

The 7-ray bursts may be accompanied by a flux of energetic neutrinos [^2|, Taking 



the distance of the 7-ray burster from the Earth to be of order 10^^ cm, one finds that all 
heavy neutrinos will decay if the lifetime of neutrino at rest, t^, satisfies the inequality 

^ lO^Sec ^. (88) 
E 

Let us evaluate this bound both for hierarchical and quasi-degenerate neutrino mass spec- 
trum setting E 1 TeV. In the hierarchical case mi ~ 0, ~ 0.05 eV and m2 ~ 0.007 
eV so from ( ^Sf ) we find that in order to let z/3 decay should be ~ 10^ sec, while for z/2, 
the bound is T2 < 10 sec. For quasi-degenerate spectrum with mi ~ m2 ~ ms = 1 eV, the 
bound is weaker: 10^ sec. 

It was estimated that the flux of neutrinos with the TeV scale energies from an 



individual gamma burster at cosmological distance ^ ~ 1 produces 10~^ — 10 muons in 1 
km^-size detectors. Since these neutrinos are correlated in time with the 7-ray bursts and 
aim at the same source, they can be distinguished from neutrinos produced by other sources. 
The rate of 7-ray bursts detectable on the Earth is ~ 10^/year so the statistics are fairly 
high and we can deduce results based on studies of such neutrinos. 

The large scale detectors cannot identify the charge of produced leptons, so in practice 
z/i and z/i signals will be summed up. Unfortunately, with present design for 1 km^-size 
detectors, it is hardly possible to identify the flavor of the detected neutrinos, for the energy 
range which we are interested. However, there are methods which open a possibility to build 
a large detector with flavor identiflcation [p^ . 

Let us assume that these technical problems will be solved and that future detectors will 
be able to discriminate between different flavors. In the presence of the decay only the flux 
of the lightest neutrino, ui, will arrive at the Earth. Then the ratio of /i-like events to r-like 
events produced by this flux equals 

/i-like events _ |f/^ip 



r— like events If^riP 
where we have taken into account that for high energies, neutrinos of different flavors have 
nearly equal cross-sections: cr(z/^) ~ cr(z/T-) and o"(t'^) — cr(z/^). 

Thus, if the detectors are able to identify r-like events, we will be able to measure the 
ratio \U^i/Uri\. Using this ratio, the unitarity condition |f/eip-|-|f/^ip-|-|f/T-iP = 1, and |f/eiP 
determined by KamLAND , one can derive the value of | f/^i | . Then | P = 1 — | t^^i P — | Ufj,^ \ ^ . 
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Similarly, for the ratio of /i— like to e— like events we have 

/i — like event s | t^^ti P 
e— like events If^eiP' 



(90) 



where we have used o"(z/^) ~ o"(z/e) and cr(z/^) ~ cr(z/e). This ratio can be used to determine 
|f/^ip immediately, once |f/eiP is known. Unfortunately, identification of e-like events is very 
difficult. 

Let us emphasize that the analysis based on (pOf ) and ( pO| ) does not depend on astro- 
physical details (neutrino production mechanism, etc.). However, one should make sure that 
all heavy neutrinos have decayed on their way to the Earth. A check can be based on ratio 
of fluxes (eventually numbers of events). If neutrinos are stable we expect Five) : Fiy^^) : 
F(z/^) ~ 1 : 1 : 1 II, while in the case of decay F{v^) : F{v^) : F(i/^) = \Uei\'^ : jt/^ip : 
If/riP ~ 1 : I : |. The above analysis was based on assumption that neutrinos from all 
7-ray bursters decay before reaching the Earth. It may happen, however, that due to spatial 
distribution of sources, the degree of decay can be different for different sources. From a 
single burst only few neutrinos can be detected, so studying neutrinos associated with only 
one 7-burst event, it is impossible to establish the existence of the decay. This can be done 
on the basis of observations of many bursts. The sources can be divided into two groups: 
close sources and far sources (the distance of the source can be measured by its redshift). 
Studying the flavor composition of neutrino fluxes from these two groups, one can check the 
stability of neutrinos. There are other measurements which can shed light on the decay rate 
of neutrinos BTI, HSi. 



Note that this analysis does not depend on the solution of the solar problem. 



4.4 Loss of coherence; averaged oscillations 

Let us consider stable or meta-stable neutrinos produced by cosmological sources. For ex- 



ample, consider again the neutrinos with ~ 1 TeV accompanying the 7-ray bursts 
For such neutrinos the oscillation length is much smaller than the distance from the source, 
L ~ 10^^ cm (even for the VAC solution of the solar neutrino problem, Arn^^j^L/E ^ 1). 
Consequently, all the oscillatory terms in the probabilities will be averaged out. Furthermore, 
according to existing models of the bursters, the neutrinos are produced in the envelope of 
the star with density p ~ 10^'' g cm~^ and radius ~ 10^^ cm and therefore the matter ef- 
fects inside the source are neghgible |H6[. As a consequence, the oscillation probabilities for 
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neutrinos — z/^) and antineutrinos {ua — > z^/?) take the form 

PaP = -Pa/3 = 51 If^mn^/Ji^- (91) 

(Here we used \Uai\ = \Uai\-) In particular, 

P'mm = E I^A'^I' = - '2K2\^\U^.l\^ (92) 



and 



= E If^/.d'|f^e.|' = " K2\\\Uel\' - \Ue2\'), (93) 



where K^^ and i^'e^^ are known functions of \Uei\, \Ue2\, \Ue3\, l^^isl and do not depend on 
and |t/^2p- The probabihty Pee does not depend on \U^i\^ and |t/^2p- 
The probabihties (pTD, (|92D, (p^) have the following properties which play a key role in 
our calculations: 

(i) -Pq/3 = Pf3a', 

(ii) probabilities for neutrinos and antineutrinos are equal; 

(iii) the probabilities do not depend on energy. 

Let us calculate the number of charged current events produced by neutrinos from 7-ray 
bursters in the detectors. We assume that the source produces (differential) fluxes of elec- 
tron neutrinos, F^, and antineutrinos, F^, as well as muon neutrinos, and antineutrinos, 
-F°, whereas the fluxes of r-neutrinos and r-antineutrinos are negligible. Using the proper- 
ties of the oscillation probabilities listed above and summing up neutrino and antineutrino 
contributions we can write for the number of yU-like events 

(/i-like events) = P'ma^ (/ P^^^dE + J F^^adE^ + ^(^^ {J ^"^ciS + / F^adE^ , (94) 

where a = (y{Uf,) ~ cr(z/^) ~ (jivr) and a = (y{Uf,) ~ o-ij)^) ~ cr(z?r) are the cross-sections 
for neutrinos and antineutrinos. Similar expressions can be written for the number of e-like 
and r-like events. Notice that the oscillation probabilities factorize out of the integrals over 
energy. 

For the ratios of event numbers we can write 

/x-like events _ P^^A + P^g 
e— like events P^e + Pe^i^ 

and 

r— like events Pru^ + Pre / x 
= -Lit — . — 1£ ^ (96) 

e— like events P^e + Pe^l^ ' 
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where 

_ / Fj^adE + J Fj^adE 
^ " / F^adE + / F^adE 

and the probabihties are defined in (|9ll) . The ratios in (|95|) and (|9^) are functions of A and 
|f/^2|- Presumably all other mixing parameters will be measured by terrestrial experiments 
described in sect. 3.1 - 3.4. The astrophysical information (and uncertainties) is contained 
in A and it will be probably difficult (if possible) to predict this quantity. So basically we 
should consider A as an unknown parameter. If future detectors are able to identify fiavor 
38| , we can determine two ratios ( p5| ) and (^) and A and \U^2\- Additional cross checks of 



results can be done if the detector is able to identify the charge of the produced lepton. 
Note that this method works for all solutions of the solar neutrino problem. 



5 Discussions and conclusions 

Reconstruction of the unitarity triangle is the way to establish CP-violation alternative to 
the one based on the direct measurements of the CP- or T- asymmetries. 

Properties of the leptonic unitarity triangles have been studied. Our estimates show that 
for maximal allowed value of jf/esl and maximal CP violation {6d = 90°) a precision better 
than 10% in measurements of the sides of the triangle will allow us to establish CP-violation 
at the 3(7 level. 

We have considered the possibility to reconstruct the triangle in future oscillation exper- 
iments. For this in the three neutrino context, one needs to measure the absolute values 
of the mixing matrix elements: \Ue2\, {Uesl, 1^^*21) l^^sl- The elements of the first side 
can be obtained from normalization. In general, the oscillation probabilities depend both 
on these absolute values and on the unknown relative phases. We suggest configurations 
of experiments: channels of neutrino oscillations, neutrino energies, baselines, and averag- 
ing conditions, for which corrections to the probabilities which depend on unknown phases 
are sufficiently small. We estimate that for value of jf/esl at the present upper bound and 
6d = 90°, the elements \Ue2\, l^esl, \U^2\, \Ufj,3\ should be measured with better than 3 - 5 % 
accuracy to establish CP-violation. 

We have found that 

1). The third side of the triangle, \U*^U^3\, can be directly measured in f^ — fe oscillations 
driven by Am^^^. In this case, the relative corrections to the probability which depend on 
the unknown phase, APs/P, are as large as |[/e3|~^e ~ 10%. This substantially restricts 
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the application of the method. Relative corrections could be suppressed at high energies: 
E > E^^ by an additional factor E^^/E. An alternative way to determine \U*^Ufj,:i\ would be 
independent measurements of \Ue3\ and |f/^3| in different experiments. 

2) . The element |f/e3| can be measured in studies of the (z/g Ve) survival probability 
at low energy (reactor) experiments. The uncontrolled corrections to the probability are of 
order of e. 

3) . I^T^sl can be determined by measurements of the v^j, — survival probability in the 
oscillations driven by Am?^^^. Both in the low energy [E < 500 MeV) and in the high 

energy {E > 5 GeV) regimes the relative uncontrolled corrections to the probability are of 
the order of e. 

4) . The ratio of elements |?7ei| and \Ue2\ can be measured rather precisely by KamLAND 
and by further solar neutrino studies. Then \Uei\ and \Ue2\ can be obtained using the 
normalization condition and the value of {Uesl (measured in other experiments). 

5) . The determination of | U^i \ and 1 11^2 \ is the most difficult part of the program. These 
quantities could be measured studying the disappearance at low energies {E < 500 MeV) 
in experiments with a base-hne L > 2000 km. The uncontrolled corrections, here, are 
relatively small (< e). 

The reconstruction of the unitarity triangle requires a series of measurements which differ 
from direct measurements of CP- and T-asymmetries. Indeed, 

• information on the absolute values of matrix elements follows mainly from studies of 
the survival probabilities; 

• both neutrino and antineutrino beams give the similar results, so that one can work 
with a neutrino beam or an antineutrino beam or with some combination of them; 

• averaging of the oscillatory terms and the loss of coherence help for determination of 
the relevant parameters. 

• the method works better (in a sense that uncontrolled corrections are smaller) for 
smaller value of ratio e = Am^„„/ Am^j^ which determines relative uncontrolled cor- 
rections. 

These factors inhibit direct observations of CP-violation. In this sense, the method of 
the unitarity triangle is complimentary to the direct determination of CP-violation from 
measurements of asymmetries. 



31 



Clearly, the proof of feasibility of the method requires further studies, and especially, 
calculations of the effects in specific experiments. In any case, realization of the program 
will not be easy. 

The ideal direct way to measure the absolute values of matrix elements would be exper- 
iments with beams of pure mass eigenstates of neutrinos. Charged current interactions of 
these beams will help to determine immediately the value of \Uai\'^ (or the ratios \Uai\'^ /\Uj3i\'^ 
if the absolute value of the neutrino flux is not known). However the masses of neutrinos are 
extremely small (much smaller than the quark case), and therefore neutrinos are produced 
and propagate in the flavor states, i.e., the coherent states of mass eigenstates. Separation 
of the mass eigenstates is a non-trivial problem. We have made some suggestions based on 
spread of the neutrino wave packets, neutrino decay, adiabatic conversion and loss of the 
coherence due to averaging of oscillations. 

In this connection we have considered several possibilities related to high energy extra- 
terrestrial neutrinos: 

1) . Neutrinos from possible annihilation of WIMP's in the center of the sun or the earth. 

2) . Decaying neutrinos from cosmological sources. 

The oscillation probabilities for stable neutrinos from far sources (for which Am^y^L/2£' » 
1) dependent only on the absolute values of the mixing matrix elements. We have considered 
the neutrinos accompanying gamma-ray bursts, as an example. 

These possibilities require large (~ Ikm^-size) underwater (ice) detectors which enable 
us to identify the flavor of the interacting neutrino. 
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Figure 1: The unitarity triangles for different values of the CP-violating phase Sr,. For 
mixing angles, we take tan^ ^12 = 0.3, sin^2^23 = 1 and sin^2^i3 = 0.12. The arcs show the 
10 % uncertainties in determination of y and z. 
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